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1. Introduction and main results 

We outline here a rather general strategy ^ to compute weak matrix elements -and more gen- 
erally operator matrix elements among hadronic states- using maximally twisted Wilson fermions 
in such a way that wrong-chirality and -parity mixings, as well as 0(a) discretization errors, are 
avoided and positivity of the fermionic determinant is guaranteed 1 . Our strategy relies on 

• the existence of renormalizable Euclidean lattice gauge models including "sea" quarks as 
well as possibly replicated "valence" quarks (and corresponding ghost fields), from which 
one can extract, among others, operator matrix elements that in the continuum limit coincide 
with those of the effective weak Hamiltonian in QCD; 

• the possibility of choosing the (sign of the) Wilson parameters for the various replica of 
(maximally twisted, single flavour) "valence" quark fields, as well as the way they enter the 
relevant quark operators, so as to avoid wrong-chirality or -parity mixings. 

The latter property has so far [jl|] been established only for the CP-conserving AS = 1,2 effective 
weak Hamiltonian operators (see discussion below) but is likely to be more general. We also 
assume in the following that the "infrared" problems in the lattice computation of matrix elements 
for multiparticle states [^] have been solved in some way, e.g. by the methods of Refs. 

We focus here on matrix elements of four quark operators. The evaluation of matrix elements 
of two quark operators without unnecessary mixings and 0(a) cutoff effects is in several cases 
feasible within the framework of plain twisted mass Lattice QCD (tmLQCD) []|] (even in the case 
of mass non-degenerate quark pairs ^) and will not be discussed here 2 . For recent reviews on 
tmLQCD see Refs. [Jsj] . An extension of our approach to the computation of Bg in the static heavy 
quark limit has been proposed in Ref. 

Other methods for the evaluation of specific four fermion matrix elements using twisted Wil- 
son fermions with no (or reduced) wrong-chirality and -parity mixings are possible and have in fact 



been developed for the cases of Bk |0|, the AI = 1/2 rule | ]TT[ ] and Bg in the static heavy quark 



limit [|12|]. For recent reviews on weak matrix elements in LQCD see e.g. Refs. []13[]. 



2. 4sNv Euclidean lattice models 

Let us start by specifying the lattice formulation of a generic (Euclidean SU(3)) gauge model 
with "sea" quarks as well as N flavours of "valence" quarks and corresponding ghosts we alluded 
to above. For definiteness we specialize to the realistic case of four non-degenerate physical quark 
flavours (u, d, s and c). The lattice action takes thus the form 



N 



S 4sNv = S g [U] + Sg/*> [Yt, %,U] +S& ,2) [Wh,Yh,U] + £ Sgf > [q f ,q f , U] + sg /2) [<j> f , U] 



(2.1) 



'We do not discuss here the various possible ways of performing the (non-perturbative) renormalization of the 
matrix elements of interest. However, if one wishes to employ a Schrodinger functional scheme without reintroducing 
O(a) bulk discretization errors, modifications of the usual setup are needed: see Refs. |Q| for two different proposals. 

2 We only recall a recently proposed application of tmLQCD to the evaluation of V us 
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where 3 S g [U] is a suitable pure gauge action , 

Stm [¥p,¥ P ,U]=a 4 Y,¥p(x) Y-V-iY 5 T:iW CT (r p )+m p -e p T 3 yr p (x), p = £,h 

X 

is the action for a pair of "sea" quarks (y/^ = (u,d), Xjfh = ( s ,c) and typically \r p \ = 1) with 

W cr (r) 



while 



5 0S 2) [9f,9f, U\ = a 4 £?/(*) r-V-iY5W C r(rf)+m f q f (x) , 

a: 

dS! 2) [*f,U\ = « 4 £0 / t Wsign(m / )r 7 .V-/ 75 W cr (r / )+m / j0 / (x) 



(2.2) 

(2-3) 

(2.4) 
(2.5) 



represent the action for a single flavour, Osterwalder-Seiler (OS) fermion [15] and that for the 
corresponding ghost field, </>/, respectively. The latter is a c-number, spin 1/2 lattice field belonging 
to the fundamental representation of the gauge group SU(3). Integration over 0y is convergent 
and yields an inverse single-flavour determinant that cancels the corresponding determinant arising 
from "integration" over the Grassmann fields qf and q~f. 

Renormalizability of the lattice 4sNv model with action ( |Q1 ) follows from standard power 
counting and symmetry [jj], [T^] arguments: for more details, such as the precise definition of the 
critical masses appearing in eqs. {22) to (2^), see Ref. [p. Here we only recall that valence flavour 
is obviously conserved and all valence (current) quark masses renormalize multiplicatively 



rhf = Z m (r f )m f , 
while the renormalized sea (current) quark masses take the form 



(2.6) 



in. 



m. 



Z p (r e )m e -Z s (r e )e e , m dxi =m 
Zp l (r h )m h - Zg\r h )e h , m Csea = m] 



-Z P 1 (r e )m e + Z s \r e )e e , 
Z p \r h )m h +Z s \r h )e h . 



(2.7) 



One can show that in a mass independent scheme all these (sea and valence) renormalized quark 
masses have a common scale dependence. 

In the following we will restrict attention to vacuum expectation values of multilocal operators 
made out exclusively of valence quark fields. Automatic O(fl0 improvement (and absence of all 
discretization errors coming with odd powers of the lattice spacing) of parity-even correlators and 
derived quantities follows from the invariance of the lattice 4sNv model under 3? x 3>j x (M — > 
—M) , with M an array made out of the 4 + N bare quark mass parameters. The precise definitions of 
3* (the physical parity operation) and ^ can be found in Ref. Since the effect of S>d x (M — > 
—M) on a generic operator of naive dimension d is to multiply it by (— l) d , besides reflecting all 
its spacetime arguments, a simple argument for automatic 0(a) improvement can be made that is 
closely analogous to the argument given in Ref. Jl7|] for plain tmLQCD 5 . 



3 For undefined notations here and in the following we refer the reader to Ref. 

4 In unquenched computations a clever choice of S g [U] can be important for the phase structure of the lattice model 
and thus crucial to be able to perform simulations at maximal twist with no metastability problems (see Refs. [|l4j 0]). 
5 A slightly different argument, which employs the r-parity properties of the critical masses, was given in Ref. rfn]. 



342/3 



Four quark operators in maximally twisted Wilson LQCD 



Roberto Frezzotti 



3. Bk and AI = 1/2 rule: prescription for bare matrix elements 

The general strategy of our approach is as follows. Let (j3 \&\oc) be the matrix element of the 
local operator G between hadronic states \a) and |/$) -evaluated in QCD with four flavours (qcd4)- 
we are interested in. Usually (j8|^|a) is extracted from a suitable qcd4 correlator, but, as argued 
below, it can equally well be extracted, up to discretization errors, from correlators of the form 



c *$l a = (*P (*uxi, ■ • •) ^(0) ®a (yi,yi,...)) 



(4sNv) 



(3.1) 



evaluated in an appropriate 4sNv lattice model where J2, <& a and <&p are suitable operators whose 
form is determined by 6 and the states \a) and |/$), respectively. Our approach is illustrated in |jl|], 
where we discuss in some detail a few cases of phenomenological relevance. The main points of 
that discussion are summarized below. One can check that in all these cases the number of valence 
(Mtm- Wilson) quark propagators to be computed is only 1.5 times larger than usually. 



3.1 K -K° mixing amplitude 

The target qcd4 matrix element (renormalized at the scale namely {K°\S'yy^ AA (iJ.)\K ) 



^jM\F^B a:(jU) , is usually extracted from the correlator (xq > 0, vo < 0) 



cXW) = ((^)W^a(o)(^)W> 



A5=2 



(3.2) 



with the bare expression of ^vv+aa § iven b y Kv+aa = (^^(sy^d) + {sy^d^sy^d) . We 
propose to consider a 4s6v model, with valence quarks u, d, d' , s, s' , c, and extract the bare matrix 
element 



16 



bare 



(3.3) 



from a correlator of the form (3.1) with <I>^ = d y$& ' , 3> a — dy$s, 



AS=2 
VV+AA 



{sy^d){s y^d') + {sy^d){s y^d') + {sy^d'){s y^d) + (sy^y 5 d')(sy^y 5 d) 



(3.4) 



and the lattice regularization fixed by 6 r c i = = r s = —r s < . Direct symmetry arguments thus 
imply |jl|] that Sv\^+AA' e 1- (^)> renormalizes multiplicatively. 



3.2 K ^ nn amplitudes 

The amplitudes for weak K — > %% decays in the AI = 1/2,3/2 channels can be expressed 
in terms of Wilson coefficients and the (renormalized) matrix elements (n + K~\d'^ A+AV (n)\K ), 
(n°n (> \^ A+AV (iJ.)\K Q ) with the bare expression of ^ y ^ +Ay (/i) given by 



'VA+AV 



(sy^u)(uy^y 5 d) ± (%d)(«y M y 5 w) 



u^> c 



+ 



+ 2L 



(sy^y 5 u)(uynd) ± (sy^y 5 d)(uy^u) 



1 








u <- 


c 


~2 







(3.5) 



The u and c valence flavours do not enter explicitly here, their regularization can thus be left unspecified. 
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We propose to consider a 4sl0v model, with valence quarks and 
extract the bare matrix elements 7 (k + K~ \j2y A+AV \K°) and (^n^^JSy^^ylK ) from (connected) 



correlators of the form ( p.l[ ), namely 

,(%10v) 
,(4sl0v) 



(3.6) 



The meson interpolating fields *!?go, ^n+n an d ^it°n are made out of valence u, d and s quarks 



only and the four fermion operator takes the form (see eq. (3.5)) 



^VA+AV 



f ±[0] 
VA+AV 



+ 



VA+AV ' 



2 VA+AV 



-ff ±[3] 
2 V VA+AV ' 



(3.8) 



with Gyf+AV involving u and c valence quark fields, 0y A }_ A y involving u' and d and so on. The key 
point to avoid wrong-chirality or -parity mixings of -2y A+AV is to fix the lattice regularization by 



rd 



-r u > 



-Tvt" 



-r r < 



r r » 



-Trill ■ 



(3.9) 



3.3 K — > 7i amplitudes 

Since in the chiral limit K — ► Tin amplitudes can be related to K — > % and K — ► vacuum am- 
plitudes, valuable information on the AI = 1/2 rule can also be obtained by evaluating in qcd4 
-for several "small" external momenta and quark mass values- the (simpler) matrix elements 
(jc + {p)\ty V+AA \K + (q)) and {p)\&y V +AA \K° {q)) , with the bare four quark operator given by 



VV+AA 



l r 
2 



{sY^u)(uy^d) ± (sY^d)(uy^u 
(sYnY5u)(uYnY5d) ± (sYnY5d){uYnY5U 



u^> c 



u^+ c 



(3.10) 



In analogy with the case of sect. (32), one can extract the bare matrix elements (n + n | ^^ +A y |^°) 



and (K°K \£±y\K°) in a 4sl0v model from (connected) correlators of the form (|3jJ), namely 



(4il0v) 
conn 



(3-11) 



where K and % is either K° and n° or K + and K + , the interpolating K and % fields involve only the 

_/9±[0] , /9±['] 1/?±P] 



m, <i and s valence quarks and, in analogy to eq. (3.8), J2, 



-VV+AA 



VV+AA 



VV+AA 2^ VV+AA 



5 &vv+aa ■ T° nave no unwanted mixings the lattice regularization can now be taken such that 



rd 



-r u > 



r,,ii 



-r r < 



-re"' . 



(3.12) 



4. Contact with QCD correlators and matrix elements 

Taking the example of the K°-K° mixing amplitude, we briefly explain the relation between 
the correlators -and derived matrix elements- in the chosen 4s6v model and those of qcd4 8 . 



7 As argued below, the renormalized counterparts of these matrix elements yield lattice estimates of the correspond- 
ingly renormalized qcd4 target matrix elements, {n + n~\ff VA+AV (jJ.)\K°) and (ifi K Q \ff VA+AV {il)\K Q ) . 

8 For more details and the analogous arguments relevant for the amplitudes discussed in sect. 3.2-3.3 see Ref. [|jp. 
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The crucial point is that the 4s6v correlator (3.1), with <t>p = d'y^s' , <t> a = dy^s and the operator 
^vv+aa specified in eq. (3.4), coincides with the qcd4 correlator (3.2), if (all) the quarks of the 
same physical flavour are regularized in the same way in the two theories. The equality immediately 
follows from the identity of Wick contractions. In particular it implies that "corresponding" local 
operators, such as ■^vv'+aa m 4s6v an< ^ ^vv+aa m <l c d4, admit the same renormalization constant 
(and thus equal non-perturbative anomalous dimension). Moreover renormalized (current) quark 
masses of the same flavour are taken equal in 4s6v and qcd4. Under these conditions, the equality 
of Cj^g4, a (x,y) and C^^^ (x,_y) carries over to the continuum limit, implying equality (as a — ► 0) 
of the renormalized matrix elements that can be extracted from them, including Bk- 

The advantage of working with the 4s6v model rather than with qcd4 is the greater flexibility in 
the choice of its lattice regularization, which results from the possibility of giving different values 
to the Wilson parameters of different replica of the same physical flavour. This flexibility has been 
exploited in sect. [O by taking = = r s = —r s i and is the key to obtain the discrete symmetries 



that ensure [JTJ] multiplicative renormalizability {i.e. no mixing) of the operator ■2w+aa- 
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